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Abstract. The spectrometers onboard the Infrared Tele- 
scope in Space (IRTS) reveal water vapor absorption in 
early M-type stars, as early as M2. Previous observations 
detected H2O vapor absorption only in stars later than 
M6, with the exception of the recent detection of H2O in 
(3 Peg (M2.5 II-III). In our sample of 108 stars, 67 stars 
have spectral types earlier than M6. The spectral types 
are established by means of their near-infrared colors on a 
statistical basis. Among the 67 stars of spectral types ear- 
lier than M6, we hnd water vapor absorption in six stars. 
The observed absorption features are interpreted using a 
local thermodynamic equilibrium model. The features are 
reasonably fitted by model spectra with excitation tem- 
peratures of 1000-1500 K and water column densities of 
5 x 10 19 to 1 x 10 20 cm -2 . These numbers imply that the 
H2O molecules are present in a region of the atmosphere, 
located above the photosphere. Furthermore, our analy- 
sis shows a good correlation between the H2O absorption 
band strength, and the mid-infrared excess due to the cir- 
cumstellar dust. We discuss the relation between the outer 
atmosphere and the mass loss. 

Key words: infrared: stars - stars: late-type - stars: at- 
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1. Introduction 

Water is one of the most abundant molecules in the atmo- 
sphere of late M-giants, and it is a dominant absorber in 
the near-infrared (near-IR) region. Water vapor in stellar 
atmospheres has been studied by theoretical and obser- 
vational methods. The strength of the H2O absorption 
possibly correlates with the spectral type, the effective 
temperature of the star, and the near-IR color (Klein- 
mann & Hall 1986; Lancon & Rocca-Volmerange 1993), 
even though such correlations were not clearly found by 
Hyland ([1974 ) . According to Spinrad 
Hyland 



spectral type M6 or later. This was consistent with hydro- 
static models of th e atm osphere of red-giants (Tsuji 1978 



p7l 
p74| ) 



Wing ( |1969D and 
H2O could only be detected in giants with 



Scargle & Strecker 1979 ). As a further complication, most 
late M-giants are long period variables, and in general, the 
band strength of the H2O absorption features depends on 
stellar variability. For example, Mira variables show very 
deep H2O absorption, and the depths of H2O features 
in Miras vary from phase to phase (Hyland 1974 ). High- 
resolution spectroscopic obs ervati ons of the Mira variable 
R Leo by Hinkle & Barnes ( 1979Q revealed that a signifi- 
cant fraction of the H 2 molecules were in a component 
with a distinct velocity, and a cooler excitation temper- 
ature than molecules near the photosphere. They inter- 
preted this 'cool component' as an overlying layer above 
the photosphere. 

These previous studies were mostly based on ground or 
airborne observations, where the terrestrial H2O interferes 
with a detailed study of the center of the stellar water 
bands. In contrast, observations from space are ideal for 
investigations of stellar H 2 features. Using the Infrared 



Space Observatory (ISO), Tsuji et al. (1997) discovered 
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a weak H 2 absorption in the early M-type star, /3 Peg 
(M2.5II-III). They argued that the observed H2O is in a 
'warm molecular layer' above the photosphere. 

In this paper, we present the results of a study of H2O 
absorption features, using data from the Infrared Tele- 
scope in Space (IRTS, Murakami et al. 1996| and references 
therein). 

2. Sample Data 

This study is based on data from the two grating spec- 
trometers onboard the IRTS: the Near-Infrared Spectrom- 
eter (NIRS), and the Mid-Infrared Spectrometer (MIRS). 
The IRTS was launched in March 1995, and it surveyed 
about 7% of the sky with four instruments during its 
26 day mission. The NIRS covers the wavelength re- 
gion from 1.43 to 2.54 and from 2.88 to 3.98 //m in 
24 channels with a spectral resolution for point sources 
of AA = 0.10 - 0.12 fim. The MIRS covers the range 
from 4.5 to 11.7 /im in 32 channels with a resolution of 
AA = 0.23 — 0.36 /im. Both spectrometers have a rect- 
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angular entrance aperture of 8' x 8'. The total number 
of detected point sources is about 50,000 for the NIRS 
(Freund et al. |1997[ ) and about 1,000 for the MIRS (Ya- 
mamura et al. 1996J). The estimated absolute calibration 
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errors are 5% for the NIRS, and 10% for the MIRS. This 
may cause systematic errors in the colors and the H2O 
index discussed in this paper. The spectra are not color 
corrected. 

All stars used in this study were observed both by the 
NIRS and the MIRS between April 9 and 24. We only in- 
clude stars at high galactic latitudes (|6| > 10°) in this 
sample, to minimize source confusion and interstellar ex- 
tinction. Each selected star has a unique association with 
the IRAS Point Source Catalog (PSC, |198S| ), within 8' 
from the nominal NIRS position, and each associated PSC 
entry has reliable 12 and 25 /im band fluxes (FQUAL=3; 
IRAS Explanatory Supplement |l988|) . Known carbon and 



S-type stars (Stephenson 1989 , 1984 ), as well as non- 
stellar objects were discarded from the sample. The cur- 
rent sample contains 108 stars. The signal-to-noise ratio 
for all stars is larger than 5 in all NIRS bands. 

Of the 108 stars used in this study, 43 have a unique 
association in the Bright Star Catalogue (BSC, Hoffleit 
& Jaschek 1991). The distribution of spectral types for 



these 43 stars are: 1 B, 2 F, 4 G, 22 K, and 14 M-types. A 
total of 31 M-giants are found in the General Catalogue 
of Variable Stars (GCVS, Kholopov et al. |l988| ), and the 
distribution of variable types in the GCVS is as follows: 
10 Miras, 13 semi-regulars (SR, SRa or SRb, hereafter 
SR), and 8 irregulars (L or Lb, hereafter L). No M-dwarfs 
or M-supergiants were found in the BSC or the GCVS 
associations. 

3. Results 

In Fig. [y, we show the composite spectra of six represen- 
tative M-giants observed by the NIRS and the MIRS. We 
indicate the position of the molecular absorption features 
due to CO (1.4, 2.3, 4.6 fim), H 2 (1.5, 1.9, 2.7, 6.2 /im), 
and SiO (8.2 /im). The broad-band emission at 9.7 /im is 
due to silicate dust. The H2O bands at 1.9 and 2.7 /im 
are visible in two early M-giants, AK Cap (M2, Lb) and 
V Hor (M5, SRb), where no H2O in the photosphere was 
expected to be detectable. 

One could argue that the early M-type stars with H2O 
absorption had spectral types later than M6 at the time 
of observation, because of their variability. However, this 
is not the case. We can estimate the spectral types, using 
the relation between the spectral type and color (Bessell 
& Brett 1988). We use the color C2.2/1.7 instead of the 
photometric color H — K, where C2.2/1.7 is defined as: 



C2.2/1.7 = log(-F 2 .2/Fi. 7 ) 



F 2 .2 and F1.7 are the IRTS/NIRS fluxes at 2.2 and 1.7 /im 
in units of Jy, respectively. For the NIRS wavelength re- 
gion, the 2.2 and 1.7 /im bands are least affected by stellar 
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Fig. 1. The combined NIRS & MIRS spectra of M-giants 
are shown. From top to bottom: (a) HR 1667 (M2III), (b) 
HR 257 (M4III) (c) X Hor (M6-M8, SRa), (d) RR Aql 
(M6e-M9, Mira), (e) AK Cap (M2, Lb) and (f) V Hor 
(M5III, SRb). The error bars represent the noise in the 
subtracted background level. Water absorption bands at 
1.9 and 2.7 /im are seen in the two early M-giants (e) 
AK Cap and (f) V Hor, in contrast to other early M-giants 
((a) and (b)). 

H2O absorption bands. Fig. [2| shows C2.2/1.7 versus the 
spectral types from BSC and GCVS of all the known K- 
and M-giants (59 giants) in the sample. There is a clear 
increase in C2.2/1.7 toward later spectral type except for 
Miras. One M9 star (KP Lyr = ADS 11423) deviates from 
this relation, but we regard this star as M5, according to 
Abt ( |1988D . 

All stars of spectral type M6 and later are above 
C2.2/1.7 = —0.085. Thus, stars bluer than —0.085 are ex- 
pected to have spectral types earlier than M6. There are 
67 stars in the < M6 region, defined by C2.2/1.7 < —0.085. 
C2.2/1.7 for AK Cap and V Hor is -0.130 and -0.102, re- 
spectively. These numbers confirm that these 2 stars have 
spectral types earlier than M6 at the time of the IRTS 
observation, even though both stars show clear H2O ab- 
sorption. 

We now discuss the relationship between the H 2 ab- 
sorption strength and the spectral type. For this, we define 
the H2O index Ih 2 o a $ follows: 



(1) Iu 2 o = log(F cont /F 1 . 9 ), 



(2) 



where -Fcont is the continuum flux level at 1.9 /im in units 
of Jy, which is evaluated by linear interpolation between 
i*i.7 and i<2.2. -F1.9 is the observed flux (Jy) at 1.9 /im. In 
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Fig. 2. The color C2.2/1.7 is plotted against the spectral 
type derived from BSC and GCVS. Ranges of spectral 
types of some SRs and Miras are represented by dotted 
bars. The thick horizontal line indicates the boundary of 
stars with spectral types < M6. 



Fig. ||, we plot /h 2 as a function of C2.2f1.7- The domi- 
nant measurement errors for C2.2/1.7 and Ih 2 o are due to 
by the slitless spectroscopy, and they are roughly 0.01 and 
0.002 for stars < M6, respectively. 

Fig. H can be used to find candidate stars (< M6) with 
H2O absorption. Since we estimated F cont by linear inter- 
polation, Ih 2 o is not zero even in the absence of H 2 0. We 
evaluate the relation between Jh 2 o and C-x. 2/1.7 for stars 
without H2O, using a linear fit for the 67 stars in the < M6 
region, by minimizing a merit function ^" \ yi — a — bxi\ 
for n data points (xi,yi) (Press et al. 1986J ). This fit is 
robust against outliers, i.e. stars with H 2 0. The result 
is indicated as a thick line. The thin lines indicate ±2cr 
level from the fit, where a is standard deviation of /h 2 o- 
Here, we use the +2a level as a threshold to find candi- 
date stars with H 2 0, because no star appears below the 
—2(7 line. The stars above the +2a are supposed to be 
stars with H2O absorptions. There are 6 such stars in the 
region of < M6. The value of 2a is almost equal to that 
of 3 standard deviations after the outliers are excluded. 
The spectra of the stars are shown in Fig. [|. These 6 stars 
above the +2a line show clear evidence of H2O absorption 
bands at 1.9 and 2.7 fim. Only AK Cap and V Hor have 
identifications of spectral type in BSC or GCVS. On the 
basis of their C2.2/1.7 (ranging from —0.113 to —0.086), 
the other four stars are probably M-type stars, and not 
K-type stars (see Fig. 0). 
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Fig. 3. Ih 2 o is plotted as a function of near-IR color 
£2.2/1.7- Symbols represent variable types in GCVS. The 
blackbody color temperatures corresponding to C2.2/1.7 
are indicated on the upper x-axis. Left of the vertical line 
at —0.085 is the region of stars < M6. The thick horizon- 
tal line indicates /h 2 o v -S- C2.2/1.7 relation for stars with 
spectral type < M6. Six stars above the upper thin line 
(+2er) with spectral types < M6 are supposed to have H 2 
absorption. Two stars from the sample lie outside of this 
plot at (C 2 .2/i.7,/h 2 o) = (0.888, 0.301), (-0.213, 0.025). 
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Fig. 4. The six early M-type stars with clear H2O ab- 
sorption bands, which are above the +2<r line in Fig. [3[ 
From top to bottom: (A) AK Cap, (B) IRAS 21222-4155, 
(C) V Hor, (D) IRAS 21269-3711, (E) IRAS 20073-1041, 
and (F) IRAS 05124-4936. 



4. Discussion 



In Fig. H, we plot Ih 2 o versus the color defined by the 
NIRS 2.2 urn band flux and IRAS 12 /xm flux (F 12 ) as: 



C 



12/2.2 



= log(F 12 /F 2 . 2 ). 



(3) 



C\ 



2/2.2 



is a measure of the IR excess due to circumstellar 



dust, and is roughly equivalent to K — [12], which is an in- 



dicator of mass-loss rate in Miras (Whitelock et al. 1994). 
Fig. H shows boundary at C12/2.2 ~ —1.0 between early 
M-type stars with H 2 and those without H 2 (we re- 
gard all stars below the +2<r line in Fig. || as early M-type 
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In contrast, our measured values for T rx and N are in 



-1.2 



-0.8 



£.4 

-"12/2.2 



0.0 



0.4 



Fig. 5. /h 2 correlates well with Ci2/2.2- The color tem- 
peratures are indicated on the upper x-axis. The open 
rectangles denote the six stars (< M6) with H2O. The 
early M-type stars without H2O are overlapped by dia- 
monds. There is a systematic dependence of C12/2.2 on 
variable types: Miras show the deepest H2O absorption, 
as well as reddest C12/2.2, an d the SRs are second. Ls are 
scattered in a regime from no-E^O to lower end of the 
H2O absorption in Miras (see also Fig. for variable- type 
dependence) . 



stars without H2O). Furthermore, there is a clear corre- 
lation between 1h 2 o and C 12 /2.2j which implies that the 
H2O absorption is related to the circumstellar dust emis- 
sion. However, as we show below, the H2O molecules are 
not necessarily located in the circumstellar envelope. 

We estimate the excitation temperature (T GX ), and the 
column density (N) of H2O molecules in early M-type 
stars. The spectrum of a representative star AK Cap is 
normalized with respect to the spectra of two early M- 
type stars, which have the same spectral types and similar 
near-IR color (C2.2/1.7), but do not show H 2 absorption 
(Ih 2 o)j or dust excess (C 12 /2.2)- The resulting normalized 
spectrum of AK Cap is fitted by a simple plane paral- 
lel model with a uniform molecular layer assuming local 
thermodynamic equilibrium (LTE) (Fig. p). The H2O line 
list is taken from Partridge & Schwenke (1997). The tur- 
bulent velocity is assumed to be 3 km s _1 . For AK Cap, 
we obtain a reasonable fit for T cx w 1000-1500 K and 
JV«5x 10 19 cm" 2 . A similar analysis for V Hor (M5III) 
results in T cx w 1000-1500 K and N ra 1 x 10 20 cm" 2 . 
If we assume that the H2O molecules were in the circum- 
stellar envelope, the mass-loss rate obtained from these 
column densities would exceed 10 -6 Mq yr _1 (assuming 
an abundance ratio of H2O/H2 = 8 x 10~ 4 ; Barlow et al. 



1996), which is about a factor of 10-100 larger than those 
expected fo r Ls and SRs (10~ 7 -10~ 8 M Q yr" 1 , Jura & 



good agreement with results by Tsuji et al. (1997), who 
suggested that the H 2 molecules in M-type stars are lo- 
cated in the layer above the photosphere. Our numbers 
T ex are also consistent with T cx sa 1150 ± 200 by Hinkle & 



Barnes (1979) for the 'cool component' of the Mira vari- 



able R Leo, which is an overlaying layer of photosphere. 
Because the H2O molecules responsible for the near-IR 
absorption cannot be in the circumstellar shell, a nd be- 
cause our results are consi stent with Tsuji et al. ( 1997 ) 
and Hinkle & Barnes (197E), we conclude that they are in 
an outer atmosphere, i.e. the layer above the photosphere, 
but below the circumstellar envelope. 

Numerical calculations of the atmospheres of Miras 
(e.g. Bowen 1988; Bessell et al. 1996| ) show that the pul- 
sation of the star extends the stellar atmosphere. Our ob- 
servations show that such an extended region could also 
be present in some SRs and Ls. The dependence of H 2 
intensity on variable types as seen in Fig. || and g may 
result from differences in the physics of pulsation. Not all 
Ls and SRs show H2O absorption, which is not surpris- 
ing, because of their complexity (e.g. Jura & Klcinmann 



1992| ). Using the light curve of V Hor (Mattei |998|), we 



confirm that V Hor is probably an SR, although it shows 
a sudden increase of its visual magnitude before the IRTS 
observation. Unfortunately, no light curve is available for 

AK Cap. 

Hinkle & Barnes ( 1979 ) found that in Miras the 



Klcinmann 1992 ). Thus, the H2O molecules responsible 
for observed feature cannot be in the circumstellar shell. 



H2O in an outer layer is responsible for the near-IR 
absorption, and we find the same situation in some early 
M-type stars. If we assume a constant abundance ratio 
of H2O/H2, then Jh 2 o is a measure of the total column 
density in the outer layer. Furthermore, C12/2.2 is a 
measure of the amount of the hot circumstellar dust, 
if the circumstellar shells of the stars in Fig. |5| have a 
similar dust composition. Therefore, the total column 
density of the outer layer correlates with the thickness of 
the circumstellar shell. This may suggest that the outer 
layer influences the mass loss of the star. 

In conclusion, we demonstrate that H2O absorption 
can be seen in early M-type stars, and that the H2O 
molecules are located in the outer atmosphere. The ob- 
served correlation between the intensity of the H2O ab- 
sorption and the mid-infrared excess implies that the ex- 
tended atmosphere is connected to the mass loss of the 
stars. 
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Fig. 6. The observed spectral profile of H2O in AK Cap 
(M2), divided by the spectra of HR 1667 (M2; open circle) 
and HR 6306 (M2; filled circle), and normalized to 1 at 
2.2 /j,m. Both HR 1667 and HR 6306 have no H 2 absorp- 
tion lines and no dust emission. The three lines indicate 
LTE model spectra with a column density of 5 x 10 19 cm~ 2 , 
and temperatures of 500, 1000, and 1500 K. The excess at 
longer wavelength is possibly due to dust emission. 
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